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Abstract 

The phase behavior of a two-dimensional colloidal system subject to a commensurate triangular 
potential is investigated. We consider the integer number of colloids in each potential minimum 
as rigid composite objects with effective discrete degrees of freedom. It is shown that there is a 
rich variety of phases including "herring bone" and "Japanese 6 in 1" phases. The ensuing phase 
diagram and phase transitions are analyzed analytically within variational mean-field theory and 
supplemented by Monte Carlo simulations. Consequences for experiments are discussed. 
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In the last decade it has been realized that soft materials can serve as versatile model 
systems to study phenomena of condensed matter physics. In particular, two-dimensional 
(2D) systems of colloidal particles interacting with light sources are ideal to mimic the 
adsorption of atoms and molecules on atomic surfaces, vortices in superconductors with 
periodic pinning arrays, and many other related phenomena 

BBS 

. Experimental studies 

on flat substrates j^] have beautifully confirmed the existence of a two-stage melting process 
mediated by the successive unbinding of dislocations [5] and disclinations [6]. If confined 
to a ID periodic potential 2D colloidal systems show even richer behavior with re-entrant 
melting [7] and novel phases such as the locked floating solid Js[ 0| . Results from analytical 
theories have been complemented by a series of numerical simulations (see e.g. Refs. [r], 
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Recent experimental [13, 



14j and theoretical (l5|, [l6| investigations have studied the effect 



of 2D periodic potentials on the phase behavior. In particular, experiments on triangular 
light lattices |l3j with a stochiometry of three colloids per site have motivated our research. 
One observes that even at rather low potential strength the colloids tend to group as trimers 
forming an almost equilateral triangle, but with a quite significant number of defects, i.e., 
groups of two or four particles. Upon increasing the potential strength the number of de- 
fects decreases rapidly such that they become unobservable already for moderate intensities. 
Then, the important low energy excitations are the orientation of the trimers, regarded as 
rigid composite objects, with respect to the lattice direction. Due to the interaction of the 
trimers long-range orientational order is expected for sufficiently strong coupling. Alignment 
of the trimers is observed experimentally [13] as soon as the defect-free regime takes over. 
Interestingly, the same authors also observed a loss of orientational long-range order at even 
higher intensities, which was interpreted as re-entrant melting and seems to be confirmed 
by computer simulations [r] ]. 

In this Letter we derive phase diagrams for commensurate colloidal systems in high 
external fields by analytic methods supplemented by Monte-Carlo (MC) simulations. The 
key idea is to reduce the problem to the low energy degrees of freedom by considering the 
integer number of colloidal particles which gather in a single potential minimum as a rigid 
composite object. Its shape is determined by the symmetry of the lattice and the number of 
the constituting colloidal particles, and its size by the interplay of interparticle repulsion and 
external potential. Short-time orientational fluctuations close to the potential minima are 
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considered to be already averaged out leaving only a discrete set of gross orientations. We 
will mainly focus on dimers on a triangular lattice since they exhibit a rich phase diagram 
and exhibit a series of intriguing phases. At the end of the letter we will shortly report how 
the phase diagrams for trimers can be obtained in a rather straightforward manner. 

For an isolated dimer there are three equivalent orientational ground states on a triangu- 
lar lattice, denoted by a = 1, 2, 3 (see Fig. |TJ). The effective interaction between the dimers 
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FIG. 1: The model system, (a) The triangular lattice of the external field and the three orienta- 
tional states of the dimers. (b) Four energy levels for the pair-wise interactions of dimers. 

results from the screened Coulomb interaction between the constituent colloids, which is 
shorted for the esperimenta, conditions in Re, Q. Therefore, if - appropriate to 
restrict the effective dimer-dimer interaction to a nearest neighbors. A pair of dimers can 
be in 3 2 configurations corresponding to one of four energy levels, with generic ordering 
Ei < E2 < E3 < E4] see Fig. ^ The relation between these model parameters and the 
experimental control parameters like screening length and potential strength can be worked 
out in detail but is far too involved to be presented here. As a guidance there are the 
following trends. Increasing the laser potential compresses the dimers and thus by increas- 
ing the effective distances between neighboring dimers reduces their interaction energies 
and hence the values of Ei and the spacings between them. This is the basic mechanism 
underlying orientational melting. Lowering the screening length effectively corresponds to 
increasing the potential strength. 

The Hamiltonian of our model system reads 

a=1 {i,j)a 

where (i,j) a denotes the pair of nearest neighbors whose bond- vector is parallel to the 
orientation of the dimer state a. (7, is a vector representation of the dimer configuration (Tj 
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at the site v. 1 i— > a = (1, 0, 0), 2 i— > a = (0, 1, 0), and 3 i— > a = (0, 0, 1). For the energy 
matrices H Q one finds in terms of a permutation matrix 
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Since a shift of the global energy scale does not affect the phase behavior only the three 
energy differences in units of k^T constitute dimensionless parameters. Thus one would be 
led to expect a three-dimensional phase diagram. However, upon rewriting the Hamiltonian 
in "spin language" it turns out that the parameter space can actually be reduced to only 



two dimensions 



17| . We find, omitting an additive constant 



(3) 



where the new energy scales K, M are given in terms of the direct dimer interaction energies 
by K = —Ei + E 2 and M = Ei — 2E 2 + 2E 3 — E 4 . This reduction greatly simplifies the 
analysis and facilitates comparison with future experiments with the two generic control 
parameters, screening length and laser intensity 22]. 

For M = the Hamiltonian reduces to a 3-state Potts model. The critical properties for 
the 2D g-state Potts model for q > 4 can be determined exactly for ferromagnetic (FM) inter- 
actions, K > 0. The critical point is known rigorously for square, triangular, and honeycomb 



lattices 



18l |l9| . There is sufficient numerical evidence that the q > 4 solution also holds for 
q = 3 although an analytic proof is still missing. The FM-P (paramagnetic) phase transi- 
tion of the 3-state Potts model is continuous as suggested by the assumed "exact" solution 
and corroborated by MC simulations and renormalization group analysis Q]. Interestingly, 
a mean-field description gives qualitatively wrong results, since it predicts a discontinuous 
transition. For negative exchange coupling K < the low temperature phase is an antifer- 
romagnetic (AM) Potts state. The results of renormalization group, series expansions, and 
MC simulations for the transition on a tri ang ular lattice appear to be controversial as far as 
the nature of the transition is concerned 19(. 

The phase diagram of dimers in the external potential forming a triangular lattice as 
prescribed by a generalized Hamiltonian, Eq. (JHJ), has a richer topology. Generically K ^ 0, 
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FIG. 2: Zero temperature ordered structures of 2D colloidal dimers: (a) ferromagnetic (FM), (b) 
herring bone (HB), (c) Potts antiferromagnetic (AM), and (d) Japanese 6 in 1 (J6/1), structure. In 
the latter, the energy of the system does not depend on the orientation of dimers on sites denoted 
by circles, thus this structure is 3 N / 4 degenerate. 

and one distinguishes between K > and K < 0. In both cases there are two dimensionless 
parameters, the ratio of the exchange couplings M = M/\K\ and the normalized temperature 
f = k B T/\K\. 

An exact analytic solution of the whole phase diagram, except for the two Potts points, 
is difficult to obtain and one has to rely on approximative and/or numerical methods. We 
start with a mean-field (MF) analysis which allows us to determine the topology of the phase 
diagram and thereby the symmetry of the order parameter in the various low temperature 
phases. In particular, we employ a variational MF approach where the full density matrix is 
approximated as a product of single-site density matrices, p({o"j}) = IL Pi(ci)- By choosing 
site-independent density matrices the order parameter of the FM phase has been found by 
minimizing the variational free energy. Appropriate generalizations have been introduced to 
obtain AM phases and other ordered phases; details will be presented elsewhere jl7 |. 

From the MF analysis we find the equation of state for the respective order parameters, 
and in particular the phase boundaries. It turns out that the variational free energy for 
colloidal dimers in the P/FM or AM phase is identical to the corresponding result for the 
Potts model, provided one substitutes the Potts exchange energy K by the effective coupling 
K — > K e ff = K + M/3. The critical temperatures for the P-FM and P-AM phase transition 
are located at ksT c = 3K efi /(2 In 2) and ksT c = — 3-K" e fj/(41n2), respectively. 

The colloidal dimer Hamiltonian, Eq. (JHJ), allows for new intriguing structures which 
are not realized within the Potts model. We have investigated for herring bone (HB) and 
Japanese 6 in 1 (J6/1) structures - a term borrowed from weaving patterns for chainmailles 



-, see Fig. The corresponding order parameters can be obtained by (numerically) mini- 
mizing the variational free energies and the phase boundaries can be determined accordingly; 
see Fig. El Except for the HB-P transition for M < —7/3 within MF theory all phase transi- 
tions are discontinuous. The HB-P transition has already been studied in the context of N 2 
absorbates on graphite. A continuous anisotropic planar rotor model analogous to our dis- 



crete version for K = yielded similar results both within MF and MC description 20l |21|. 
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FIG. 3: Phase diagram for dimers. Solid/dashed lines represent discontinuous/continuous MF 
transitions. Dots and crosses denote the continuous and discontinuous transition points as obtained 
by MC; the two encircled symbols indicate the pure Potts transitions. The nature of the HB-P 
transition for M < —4 is not clear which is indicated by a dot-cross symbol for the critical point. 
Dotted lines represent extrapolations of the Potts model critical point, see text. 



In our case, the MF description is not expected to predict the correct order of the tran- 
sition as can be inferred by specializing to the Potts model (M = 0). Nevertheless, the 
MF results should constitute reasonable approximations for the actual phase behavior. In 
order to gain further insight we have performed extensive MC simulations of the colloidal 
dimer Hamiltonian, Eq. Q, using a standard Metropolis procedure. The observables have 
been averaged over > 5000 statistically independent configurations after the system has 
been equilibrated. Every 1000 MC cycles a new configuration contributing to the averages 
is considered. Simulations have been performed on several system sizes, N = L x L and 
L = 12, 36, 54, 108, and periodic boundary conditions have been employed to mimic bulk-like 
conditions [ljj. The order parameter characterizing the various phases has been chosen in 
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analogy with the Potts model, either defined on the whole system or suitable sublattices 
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Here, N a is the number of dimers in state a, in particular, iVs corresponds to majority ori- 
entation, and iV = J2a N a is the total number of the dimers on the whole (sub)lattice. We 
have performed temperature sweeps for K > and K < with a number of choices for M 
starting from the ordered structures. We have monitored the average order parameter and 
energy, as well as their respective variance. The transition temperatures corresponding to 
melting of the ordered structures have been determined by extrapolating the location of the 
maxima of the heat capacity and susceptibility to the infinite system size. In order to decide 
whether the respective phase transition is strongly discontinuous or compatible with a con- 
tinuous/weakly discontinuous transition scenario we have evaluated Binder's cumulant and 
respective probability distributions for the order parameter and energy both scenarios 
are exemplified in Fig. 0] 
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FIG. 4: Order parameter 5 as a function of the reduced temperature for the two FM-P transitions. 
Left: Potts-like continuous transition for K,M > 0. Right: the M-term driven first order transition 
for K < 0, M > 5. Inset: probability distributions of states with respect to the energy per dimer, 
u = U/N, below, above, and at the transition temperature. 



Figure 01 exhibits the phase diagram for the colloidal dimer Hamiltonian; boundaries 
resulting from both MF description and MC simulations are included in the Figure. The 
variational MF free energy reproduces the topology of the phase diagram, although some of 
the transitions become continuous by fluctuation effects. Quantitatively, MF overestimates 
the transition temperatures with respect to simulation results by a factor of order unity, 
a feature already familiar from the Ising model Interestingly, the discontinuous HB- 
P transition appears to saturate for large negative values of M in strong contrast to the 
corresponding MF result j^j. For colloidal systems the regime corresponding to K > 



and M < should be experimentally accessible. A stringent test of our theory would be to 
experimentally verify the predicted HB-P phase transition. 

Fluctuation effects reduce the critical temperature of the Potts model by a known fac- 
tor [19]. The variational MF free energy of the colloidal dimer Hamiltonian for P/FM and 
AM structures is equivalent to the MF Potts model, provided that the Potts exchange energy 
is replaced by the effective coupling K e $. This analogy suggests an empirical improvement 
of the MF result for these transitions by rescaling the transition temperatures with the same 
respective correction factors. The result of this procedure is shown as dashed lines in Fig. 01 
For K > there is excellent agreement with simulation results for temperatures T > 1. For 
K < the mapping gives an overall quantitative improvement of the phase diagram. In 
particular, for M = the well known Potts solution is recovered. 

We close with reporting our main results for trimers on a triangular lattice. Due to the 
symmetry of the lattice, trimers can only be in one of two orientational states, represented 
by Ising spins a = ±. Hence the four configurational states for trimer pairs contribute one 
out of the three possible interaction energies, 

mm" 

+ + 

Ei < E 2 < E 3 . (5) 

In contrast to the Ising model the interaction energies for the (H — ) and ( — h) configura- 
tions are different. The relation between the model parameters and the control parameters 
of the experiments are similar as for dimers. In particular, an increase of the potential 
strength again leads to a reduction of the interaction energies between trimers and thus 
favors orientational melting. 

The central questions now are how many relevant energies determine the phase behavior 
and what is the effective statistical model. Similar considerations as for dimers would suggest 
that there are two relevant energy scales resulting in a two-dimensional phase diagram. 
Actually, by a simple rearrangement the Hamiltonian can be cast in the form of a Ising 
model with a single excitation energy AE = E\ + E 3 — 2E 2 corresponding to single 
trimer flips, irrespective of the configurations of the neighboring trimers. Depending on the 
sign of AE the ground state is either "ferromagnetic" with parallel (++) (or equivalently 
( )) configurations or "antiferromagnetic" with an alternating arrangement of ( — h) and 
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(H — ) configurations. It is interesting to note that it is the relative orientation of trimer 
pairs and not the trimers itself which determines the ground state. 

For the remainder of the analysis known exact results of the 2D Ising model Q| can now 
be employed to obtain the phase diagram. Considering the parameters of the experimental 
set-up in Ref. [llj we can identify the critical laser intensity, V = 78 k B T Q], above 
which the orientational order of trimers is lost. This result is consistent with experimental 
observations, where the colloidal trimer system was found to be orient at ionally ordered for 
V = 60 k B T but disordered for V = 110 k B T. It would be interesting to extend the 
experimental analysis to find the actual transition point and compare it with our prediction. 

We thank C. Bechinger for discussions. A. S. acknowledges the support by the Humboldt 
Foundation. 
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